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ABSTRACT 

A theory of electron-molecule scattering is developed which is a 
synthesis of close coupling and adiabatic-nuclei theories. Specifically 
the theory is close coupling with respect to vibrational degrees of free- 
dom but adiabatic-nuclei with respect to rotation. In addition, this 
theory can be applied to any number of partial waves required, the re- 
maining ones can be calculated purely in one or the other approxima- 
tion, A theoretical criterion based on fixed-nuclei calculations and 
not on experiment can be given as to which partial waves and energy 
domains require the various approximations. The theory allows all 
cross sections (i.e., pure rotational, vibrational, simultaneous 
vibration- rotation, differential and total) to be calculated. Explicit 
formulae for all these cross sections are given. 


"“NRC-NASA Resident Research Associate. 



The theory is applied to low energy e-N, scattering. The fixed- 
nuclei results are such that the criterion shows clearly that vibrational 
close coupling is necessary, but only for the 11 partial wave. The con- 
tribution of remaining partial waves can be obtained directly from the 
adiabatic-nuclei approximation. The close coupling calculation for tlie 
rig wave is carried out, and we find that it does give rise to the sub- 
structure as well as the gross structure of the 2.4 eV resommce. When 
tliis amplitude is combined with the adiabatic amplitudes we can compute 
absolute values of all cross sections of interest. In particular we find 
that vibrational excitation cross sections ax’e about twice as large as 
previously inferred. The momentum transfer cross section can also be 
computed, and it too reveals substructure witiiin the gross structure 
resonance. 
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I. INTRODUCTION 


Molecular nitrogen is a major constituent of the atmosphere to an al- 
titude of about 500 Ion. Thus it can be expected that the scattering of elec- 
trons from N 2 will be an important process in the aeronomy of the atmos- 
phere particularly above the E region, where photo ionization of many of the 
upper atmospheric constituents Iiy solar UV produces an abundance of elec- 
trons. For example the resonances in the e-Nj vibrational excitation cross 
sections in the vicinity of 2,4 eV have been user] by Newton et al,,^ to ex- 
plain an enhianced + N 2 (v) NO"^ + N rate which in turn will lead to the ob- 
served decrease in ambient electron density in the Fj region during times of 
enhanced air-glow giving rise to stalaLe auroral rod (SAR) arcs via the reaction 
e + NO+ ^ N + O. 

The 2,4 eV, e-N 2 resonance is known to bo a veiy complex structure, 
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Fortunately there is a wealth of experimental details, and more or less 

phenomonological theories to explain them. It is important howover that this 

complicated structure be understood from a fundamental— essentially ab initio 

point of view, in order that researchers be able to predict scattering from 

other molecules which cannot be prepared in the laboratory, but which we 

3 

now Itnow to exist in many astrophysical environments. We telicve that the 
present modification and synthesis of existing theories will complete the fund- 
amental approximations which must underlie the methodologies to be employed 
in such calculations. 
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i 


In Section II we fshall describe tJie e-Nj, 2.4 eV resonance and briefly review 

the previous calculationiil theories leading to the ijoint in Section JJJ tJiat strictly on 

the basis of I'ixcd-nuclei calculations (and not exjjeriment) of the JI^ parliaJ wave, 

one ojin infer that the adiabatic-nuclei theory will not sulTice to describe the sub- 

4 

struedure of the resonance, i’he fixed-nuclei calculations which are an extension 
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of those of Burke and Chandra to a series of internuclear distances H are also 
described in Section III. In Section IV we develop the vibrational-close coupling 
theory and show that our calculations for the partial wave do reveal the sub- 
structures of this resonance as more states are added while at the same time 
showing reasonable convergence when a sufficient number of states is retained. 
Tlie main formulae of the adiabatic- nuclei approximation are reviewed in Sec- 
tion V; by simple inspection one may then see how to combine the vibrational 
close coupling with adiabatic-nuclei theories: the process is one of substitution 
of the appropriate scattering matrices of the one theory for the corresponding 
ones of the other. 

Results and comparisons with experiment are also given in Section V. An 
important aspect of the ealculatio is is that ths 3 ' yield absolute normalization for 
individual vibrational excitation cross sections, whereas as far as we know all 
measurements are relative or have been done at indivdduai angles. Furthermore 
the theory is complete: it gives formulae for total and differential cross sec- 
tions of individual and averaged transitions for vibrational and or rotational 
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excitation ns ivell as pure clastic scattering. Simultaneous NTbraticn, ivtation 
differential cross sections can also be calculated (although that is not ione herd 
as well as the momcntuiu-transfcr. The latter, wiuch is calculated, is particu- 
larly interesting because it reveals substructure similar to that of the elastic 
scattering. 


Finally in Section VI we give a brief discussion of how this method fits in 
with other methods and possible generali/.aiioits. 


n. PRELIMINARIKS 

Experimentally the low energy o~Nj resonance is a oomplicaiod beast. It 
was first measured in detail in vibrational oxeibition by Sehul?: as a series of 
irregular peaks (cf. Figs, 11, 12 & 13). The rcstilts have been imporbuitly 
complemented by the ineasurGmout of the vibraWonally elastic scaitering hy 
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Golden, which shows a gross stmeturo peak centered at about 2. teV super- 
imposed on top of at least five prominent sub poal^s between 1.8 and 3.2eV (c.f. 
Fig. 7). 'I’hc observed inelastic structure of the led in short order to its inter- 
pretation (crudely dose ribod) as a compound state ofthc electron and tai'get system 
8 9 

i.e, , the ion) ’ with the substructure due to the interference between the 
various vibrational states of the compound system. This general physical in- 
terpretation of the Lstructuro is ceiiainly correct; the difficulty with the calcula- 
tions is that they are not ab iulHo , and thus they contain many adjustable param- 
eters. It should be added, however, that the physical and mathematical refine- 


ments of this general approach have been greatly developed since the original 


papers, most successfully by Birtwistle and Herzenberg. That calculation which 

is based on a boomerang model, gives remaricable agreement with the slmpes of 

11 

the vibrational excitation curves as measured by Ehrhardt and Wilimann but it too 
does not give absolute values. 

In addition to the types of calculationB mentioned above, there have been 

other, more ab initio types of investigations. One Is a calculation by Krauss and 
12 

Mies of N' as a bound state structure. This calculation confirmed the assign- 
13 

ment by Gilmore of this resonance as a fl "state," and it showed that it is a 

a 

shape rather than a Feshbach resonance. The calculation is truly a 15 electron 
self-consistent field calculation, but since the Ilg symmetry overlaps the ordi- 
nary e-Nj continuum, which can be of lower energy, some delicate restrictions on 
the variation were necessary to assure that the resonant state was not contaminated 

by this nonresonant scattering of the same symmetry. On that score we can all 

12 

state that the calctdation was in good hands with the NBS investigators. 

Finally two single-center fixed-nuclei calculations have been carried out by 

14 

the Belfast group: the first by Burke and Sinfailam includes full exchange of 
tlie incident and orbital electrons, but no induced polarization; the second by 

g 

Burke and Chandra includes polarization and simulates exchange by orthogonal- 
izing the scattered to the bound orbitals. This 'pseudo-potential^ approach pro- 
vides the basis of all fixed-nuclei aspects of the present calculation, and we shall 
discuss it as appropriate in succeeding sections. 
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m. PKED-NTUCLEI THEORY AND CALCULATIONS 


In contrast to its application in bound state problems, the fixed- nuclei ap- 

X5 1 6 

proximation for electron-molecule scattering * assumes not only that the 
nuclei are fixed (at a distance R apart) , but tliat the target molecular wave func- 
tion tI>(x;R) has been precalculated at each, in principle arbitrary, internuclear 
separation R. As a result the scattering associated with a total wave function 
(r is the coordinate of the scattered and Xj those of the orbital electrons) 

( 3 . 1 ) 

exhibits no particular stationary properties with respect to variations of R about 
the equilibrium separation of the target molecule R = R^^ . This is a very force- 
fully exhibited here, when we extend the calculations at the equilibrium separa- 
tion. (R = Rq = 2. 068a(j) to four additional values of R of the Nj wave function 
given by Nesbet:^'^ R = 1.744393, 1.868, 2.268, 2.391607, 

To review the fixed-nuclei calculation (cf . Ref. 5 but our notation is some- 
what different), the ground state * 2g of N 2 is a closed shell described by a single 
Slater determinant 

^ = det (0ai(xi)0a2 ^^2) • • - ^>0:14 (X14)) ( 3 . 2 ) 

where the Oj can readily identified from the configuration lUg 2 ct| 30g 

17 

20 y iTT^ of the ground, ^ Sg , state of N^. The bound state function which is 
an LCAO function- meaning pairs of orbitals are centered about the separate 
nuclei~is converted to a single center basis using a program of Faisal and 
Tench. 
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WO' T, -«s“(r)Y«(n). 


(3.3) 


(A double prime indicates every second term is to be taken. About eight terms in 

this expansion are sufEcient for convergence. ) A similar expansion is now made 

IG 

of the scattered orbital 


'I'raW = E 

e 






(3.4) 


where are the sperical angles of the intemuclear axis in the laboratory frame 
and r = (r, i.e, , unprimed coordinates are the coordinates in the molecular 
frame. 


With the use of (3 . 3) the static potential seen by the scattered electron is 
naturally expanded in single center coordinates (in rydberg units) 

N 




2Z 2Z 


-E 




= E VxCrtPxCcosO) (3.5) 


|r-R/2| |r + R/2l |r-rj 

The equations satisfied are then derived from the variational principle; 

6 / * (H - E) dr = 0 (3.6a) 


winch is equivalent to the projection 

/Y?,™Yc,„(H-E).I.™'dr-' = 0 (3.6b) 

wher'^ dr“^ means integi-ation over ail coordinates but r (including integration 
over ). In practice this set of coupled equations, which may readily be 
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derived fi’CTn (3,1) - (3,5), is augmented to include an induced polarization 
potential 


{ \\ /<^n(R) Oll(R) 

V(P"l)(r;R) = - ^ 

\ *■> 


P,(COS0) [l-e-Wro)*l, 


(3.7) 


where is greater of r and VaR. The calculation also includes orthogonality 
to all occupied orbitals of the same symrketry via Lagrange multipliers. The 
equations satisfied by of Equation (3.4) are then 

ug!(r). S 4”,'(r)ugl(r) = E ^ <3-B) 
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> J 


1 J 


where 


/ 2JJ, _ I 

^ ^ (2j X o o/J2j o)(Cj X m o/fij m) Vj^(r), 


(3,9) 


and 


v;^(r) = Vj^(r) + V^P''->(r). (3.10) 

Y^pob multiple components of from (3.7), thus in particular 

0 forX>2. [The remaining sjonbols in (3.9) are Clebsch-Gordan coefficients,] 


The calculation for each R was done just as the calculations of Burke and 

C 

Chandra atR==Rjj. One only needs the dependence of the polarizabilities on R. 
These were taken of the form 


a^CR) = 12.0+ 1.692(R-Ro) 
oc^(R) = 4.2 + 2.03 KR-Rq) 


(3,1 la) 
(3.11b) 


g 


Eqiiations (3.11) were chosen to give the correct polarizabilities at R = and 

to reduce correctly to the united atom limit: Oo(0) [= Osilicon^ =8,5aQand 

19 

=0* The value 8 . SSq was interpolated from Stemheimer's calculation 

of the polarizabilities of C2", K'*', and Ca'*”^. It is to be noted tlmt Equations 

20 

(3.11) are somewhat different from Truhlar who used Raman data to get an 
accurate estimate of the derivatives in the neighborhood of R = Rq. Our own 
interpolations, while somewhat cruder, should apply over a larger range in R, 
and thus be more suitable to excitation of higher lying vibrational states, Einally 
the value of was retained at 1. 592 as independent of R, That value was 
chosen, so that the Ilg resonance for R = R^ occured at exactly k2=2.394eV. 


The scatteiing is determined from the asymptotic solution of (3.8): 


lim 

r CO 


“y* = k>4|siii(kr- 


* Kb's’ 

i j 



(3.9) 


The K matrix as indicated in (3.9) is diagonal in m; it is also real, ssnnmet- 


ric, and for homonuclear targets connects only 2j and 2^ of the same parity. In 
matrix notation the scattering is naturally expressed in terms of a matrix pro- 
portional to the T matrix. In Reference 5 tliis matrix is taken to be 
which is related to by 

T(m) = 2i(l (3.10) 


In Reference 16, which gives the original derivation of the coupled fixed nuclei 

cross sections (using a spherical analysis), the scattering is written in terms of 
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the a matrix which is related to the above T matrix by 




j 

i 




^ .Ci-Ci-1 V^r(2Gj + i) 

* '* — — Ann 


C-fi- 
I j 


(3.11) 


The fixed-nuclei results are most conveniently given in terms of the sum 
eigenphase shifts. The eigenphases are the arc tangents of the eigenvalues 
of the K matrix: 


detiK^^iJ II = 0 (3.12a) 

^ tan-i (3.12b) 

i 

where the STim in (3.12b) goes over all coupled states that are included (and is 
found to converge with inclusion of approximately 8 coupled states as stated 
above.) | is the unit matrix, in Table I we give a selection of our results for 
, Zy, partial waves as a function of intemuclear separation R, The 

g 

R(=Rq) = 2.068 results are just those of Burke and Chandra. A detailed de- 

22 

scription of that generic program has been published by one of us and that 
program is what was applied here. 


We also note in addition to m = 0,1,... corresponding to Z , IT, . , . that 

tlie parity of the index 2 (even or odd corresponding to g or u) is also a good 

quantum number as well as the spin S, The latter is always S = Vi corresponding 

to doublet multiplicity, since is a closed shell (^Zg) target. [We therefore 

suppress the doublet label, for example , on our partial wave notation.] 

& 


The sum of eigenphases for non-n phase shifts are seen to change min- 

b 

imally as a function of R (although it is interesting that for Zg the minute change is 
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an oscillatory one). The change is also slow and smooth as function of the im- 
pacting energy k^. For those partial waves therefore the adiabatic-nuclei theory 
for both rotational and vibrational excitation applies (see below). 

On the other hand the change witli both R and of the Ilg wave, given in 
Figure 1, is dramatic! As a function of k^ the subient feature is the resonant 
behaviour. If one confines attention to the equilibrium separation R , one evalu- 
ates the wldtli r = 0,4eV, wliich is sensibly larger than the vibrational spacing 
AEy = 0,29eV. It was for this reason that we previously believed the adiabatic 

nuclei theory would be at least semi-quantitatively applicable to that partial wave 
23 

as well. However if one looks at the curves for R > R^ then one sees that the 

resonance has diminished to F = 0.14eV for R = 2.391607. And even at R=2.268, 

r = 0.25 eV which is smaller than the vibrational spacing of Nj; in other words 

the time (r cc ) spent by the incoming electron in the vicinity of the molecule 

is comparable to or longer than the vibrational period of the nuclei. TMs is a 

definite violation of a basic criterion for the validity of the adiabatic-nuclei 

theory, and it gives a purely theoretically determined basis for distrusting the 

24- 

adiabatic-nuclei theory for tins paiiial wave. We therefore turn in the next 
section to vibrational close coupling and its amalgamation into the adiabatic- 
nuclei theory. 

Before concluding tills section, we give in Figure 2 a comparison of our 
fixed-nuclei width and position curves vs. R as compared to the calculations of 
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Krauss and Mies and Birtwisth and Herzenberg, Considering the different 

natures of these calculations, we consider tlie agreement to be remarkable. 


IV. VIBRATIONAL CLOSE COtFPLING 


It is clear from the foregoing that it is necessary to include the dynamical 
response of the nuclei to their vibrational motion. The most natural way of do- 
ing that in quantum mechanics is to expand the wave fjnetion in terms of the 
eigenfunctions of the vibratioml motion: this is what is meant by a vibrational 
close coupling expansion; 


•I, cm) = ^ X^CR) (4.1) 

V 

[The contrast of this with (3,1) should be noted.] Let us write the total 
Hamiltonian (in rydbergs) 


H = H„{R) - - 1 + V,, + V„ , (4.2) 

where Hq(R) is the Hamiltonian the target molecule, with nuclei at a distance R 
apart and M their reduced mass. Vg^ and V^e are the interaction potentials of 
the scattered electron with the or'bital electrons and nuclei respectively: 

N 


Y = V' 

ee ^ _ r \ 


(4.3) 


1 = 1 


Vn, = -2Z 


(4.4) 


^lr-R/21 1r + R/2iy 

We derive coupled equations for the functions F^(r) in the usual way; obtaining: 
[-VJ -kj,] F„,(r) + 2] <®''' + F,(r) = 0 (4.5) 
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In. deriving (4. 5} one uses the conservation of energy 


E-E„ = k2 = k;+e^ (4.6) 

where is the electronic energy of the target state satisi^lng the target 
Schro dinger equation; 


H^CR) R) = EqCR) ^h(rj, R) (4.7) 

One also uses the approsdmation that the rotational Idnetic energy is negligible 


compared to its vibrational energy, so that 

1 d^R 
^ 

^ R dR^ 


(4.8) 


and the fact 


[- 


1 1 d^R 


+ (ch(rj, .R) iH^(R)|<I>(r;, R)\-e^ 


Xv(R) = 0 


(4.9) 


M R dR^ -V 

In (4. 5) and (4. 9) the subscripts on < > indicate the coordinates over which one 

integrates. Note in particular that <<lV ^ includes the parametric 

dependence of on R and is therefore n^ simply <4> i V^g 1 4>)j. . On the 

other hand it is true that 


®v' |V„,1 = <v' IV^^I V>K (4.10) 

since 4> is normalized for each R and V^^g is independent of I’j . The net results 
is that the set of equations, (4. 5) is seen to be a set of equations purely in r. 


To eliminate the angular dependence we malce the usual type spherical 
harmonic expansion. As opposed to (3.4), however, we here suppress the de- 
pendence on , since this does not alter the dynamical equations. Let 
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(4.11) 


= 7 E ft” 


«>IV„e+V„J<I»,. = E vW(r. R) P^Ccos «) 

‘ A 


. The latter implies that 


= 2]vf^^)(r)Px(cos0) 


(4.12) 


(4.13) 


where 


V(^^)(r) - <v'lvC^)(r,R)jv>j^ 


(4.14) 


Using all these plus the well Imown formula for the integral of three spherical 
25 


harmonics 


= (-^'"[(^S' + D(2fi' + 1)]'^ 


one can reduce (4.5) to the coupled equations: 


9! X £\/fi'Xe' 


.m o-m/ \o o o, 


(4.15) 


£! 

dr^ 


C'(C' + 1) 




= E V(28 + 1)(2!!' + 1) (-1) 

V, 2, X 


m 


'je' X A/e'xr 


vm o -m/ \ o o oy 


(4.16) 


"J'vXrt ttf’w 


Equation (4, 16) shows that the price we have to pay to get rid of the angles in 
(4,5) is the fi coupling. In addition to that we have the summation over X wliich 
comes from the single-center spherical expansions of Vg^ and Vj^g, 
Equations (4.3), (4.4). Finally however in the spirit of our fixed-nuclei 
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calculation we include an additional pseudo-potential to describe and simulate 
the effects of polarization and exchange. [Note, like (3.1), (4.1) is also not 
anti-symmetrized between r and (i = 1, 2, . . , N)1 . However the way to do 
is now clear; to (r) we add a polarization potential gotten from appropriate 
matrix elements of (r, R) of (3.7). i.e. , in (4,5) we augment the potential 

by <'hv' d’v>; tliis induces a change in the potential of (4. 14) to T^^J(r) 

in Equation (4,16), where 

*’ ?v’ ' + <''■ I ''V W- ' 

and is the multiple component of [see below Equation (3.10) J. 


Equations (4,16) are solved in analogous fasliion to (3.8) specifically here 
one demands the asymptotic form 




( jr2 

vgcosjkyr-y 


14.18) 


From the matrix one can, in analogy with (3.10) and (3.11), develop a 


and or matrix 


?'Jr,v! = 2> E 

II ,11 ’ 


vC 


( 4 . 19 ) 


V , x: 


or 


tm) _ 1) „(m) 

..^n* _.n " ^ 


vT, v2 


v'2',vS! 


(4.20) 
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The actual calculation, as we have indicated, need only be done for the Fig 
partial wave. For this partial wave there is no orthogonality requirement, on 
the other hand the multiplicily of coupling (v, 2, and X) mal^es it impossible, 
even on our macliine (IBM 360-91), to include a sufficient number of terms to 
get Ml convergence in all coupling indices. We have therefore chosen to de- 
limit the C coupling to three terms (2=2, 4, 6), The X expansion is thereby 
automatically restricted to seven terms (X = 0, 2, . . 12). Witiiin 

this approximation we seek convergence in v . The role of the polar- 
ization potential , specifically the cut-off r^ here will serve, in addition 
to exchange, to simulate the unincluded £ and X components. It was chosen so 
that with the inclusion of only one vibrational state, the resonance in the 
V = 0 V = 0 cross section occured at = 2.4eV (cf. Fig. 4). 

We shall not dwell on the numerical aspects of this calculation: suffice it to 

say that the generic program of Reference 22 was applicable with only minor 

additions to (4.16). The static potentials were generated numerically 

17 

from Ncsbet's wave hmctions; a selection of these potentials as a function of 
r is shovTi in Figure 3. Tliere diagonal potentials are compared with the fixed 
nuclei (static) potentials V^^^(r, R^) wherein we see that both potentials get 
increasingly more sharply pealced and concentrated around r = Rq/ 2, but that 
for coi'respondiug X the close coupling potentials are softer and without cusps. 

The off diagonal potentials have no real counterpart in the adiabatic-iiuclei theory, 
and they are mathematically the source of the substructure of the Hg resonance. 
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To see how this sTihstructure appears we show in Figure 4 the IT contribu- 

o 

tion to the V - 0 “>■ v' =0 cross section. The cross section is plotted for differ- 
ent numbers of vibrational states retained in the expansion. In addition we have 
shown two sets of curves, one includes two coupled 2(2 =2, 4) components in 
(4.16) and one set includes three 2’s (2, 4, 6), For each case we include all X 
allowed by vector coupling. 

One can see that one v term result does indeed exhibit a resonance at 
= 2.4eV, It was essential to get tins resontince that the polarization potential 
be included along with the static potential in , To get the resonance at the 
desix’ed position we had to chose = 1.490, 1,554 for two— and three 2-coupled 
calculations respectively. These values are gratifyingly close to the value needed 

g 

in the fixed-nuclei cEilculation 1.597. Thereafter one sees as the fundamental 
result of this paper that the substructure begins to appear, and that by the time 
we have coupled in 10 and 9 vibrational states respectively, reasonable (but not 
precision) convergence is seen to occur. (The program correctly includes 
whether various vibrational channels are energetically closed or open.) 

The comparison of the two vs, three coupled 2 solutions shows that detailed ef- 
fects do depend on the number of 2^3 retained. Thus to obtain the second rather 
than the first bump as dominant in as is revealed by experiment (of. Fig, 13), 
the retention three 2’s is necessary. It is also clear that the inclusion of even 
more 2’s would be necessary to obtain the totality of peaks in the substructure. 
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(We shall see other mauifestations of the truncated substructure in other cross 
sectioned data as well.) 

In Figures 5 and 6 we give similar results for the (Hg contribution) to 
V = 0 ^ 1 and V = 0 ^ 2 cross sections. The same type of convergence and 2- 
conpling effects are apparent, except that the necessity of more £- coupling (and 
consequently more v coupling) becomes progressively greater, as one goes to 
higher v'. 

We also state for the i^ecord that the substructure does not appear if we re- 

23 

tain only one 2. Nor does it occur in the adiabatic-nuclei approximation (see 
below) . 

V. ADIABATIC-NUCLEI APPROXIMATION, HYBRID THEORY, RESULTS 

26 

The adiabatic-nuclei approximation gives the transition amplitude between 

vibrational, rotational states P P' as 

<r'if(|,n')|P>5+e (5.1) 

where f(|, is the fixed- nuclei amplitude with nuclear coordinates frozen 

at Although the error term e on the RHS of (5.1) has never been completely 

elucidated, the present calculation will show that the delay time of scattering 

must be small compared to the period TShj/AEp'p associated with the largest 

energy quantum number of P which changes in the transition (we assume P’ # P), 

in order for e to be neglibible. Hie fixed-nuclei amplitude for scattering from 
0 

a diatomic molecule is such the angular integrations in (4.1) can be done an- 
27 

alytically, whereas the integrations over R are necessarily numerical. 
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J 


The resulting formula for the simultaneous rotation-vibration differential 


cross section (when the target is a Sg state) can be written 


28 


dajVjv _ kj'v' I 


da' 


'JV 


s A^CiV; jv) Pl(cos 6') 
47T l 


(5.2a> 


where 


AlG'v'; jv) = (2L + l)(2j' +1)2 7(2J2i-M)(2Xj-i-l) 

\o o o/\o o 0/ J \m-m o/\ JU -/i o / 


/j' j jyre. X, p 

\o o o/ |Xj Cj 
29 

and for joint vibration- rotational excitation 


“2=2,111^^'’ V) = r Xv'(R) ag £ (R) Xv(R) R^ dR 

‘ J Jo ^ ^ 


(5.3) 


If one averages over initial rotational states j and sums over final states j', one 
arrives at 


da, 


V , V _ ^ 


^ — S A, (v',v)Pl(cos0) 


dn 47 t 


(5.4a) 


where the rotationally averaged coefficients are; 
A,^(v’, V) = (2L+ l)2u2«i + IKZXj +1)1’^ 


(5.4b) 


'S, \ LU^ Xj L\/e, Xi L 


■Sj Xj L 


P, (cos O') 


\o o 0/^0 o o/\m [1 -(m+ja)/\m fi -(m+/a)^ 

The total cross are quite obviously just 4 tt times the L = 0 terms of the respec 
tive expressions 

_ kj.,. 


^jVjv = f- ^o(jv;jv) 
kjv 


^v',v = 


ky ^(2X+ 1) 


(5.5) 


(5.6) 
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[The foiTHulas involving the sum- averaging over j states are the ones wliich 
are presently useful in comparison with experiment in e-N^ scattering, because 
the rotational spacing lias not as yet been resolved. ] 


We now come to the fundamental statement of the hybrid-theory; replace 
the abiabatlc-nuclei matrix elements . Equation (5.3), by the corresponding close 
coupling values. Equation (4,20), for whatever partial waves are necessary; 




(m) 
v'E, vX 


(5.7) 


(An equivalent replacement can be made for T and K matrices also.) The justi- 
fication for this replacement we hope is clear. The point, we wish to reemphasize 
is that the time delay criterion can be theoretically (and need not be experimentally) 
assessed. (It should be added, however, that any such criterion is always aph- 
proximate in the sense that constant of proportionali'fy is necessarily somewhat 
ambiguvius. In addition, in the present case the Rq = 2.068 curve lies slightly 
removed from the main portion of the wave function, (R), of the zero vibra- 
tional function of the N 2 molecule. Thus the theoretical nature of the time delay 
criterion must be understood within the confines of such mundane considerations.) 


Coining back to the hybrid theory, note that the purely rotational aspects of 

even this Ilg wave are described by vector coupling coefficients which derive from 

the adiabatic-nuclei theory. It is this fact which has motivated us to call this a 

hybrid theory. One could in principle also include rotational close coupling at 

30 

the same time as vibrational close coupling (cf. Discussion). Indeed Henry has 
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attempted that in the case of e-Hj scattering with only limited success. The 

28 31 33 34 

point is that the adiabatic -nuclei theory is quite sufficient* ’ ’ ’ when 

polarization is included, to explain the e:iq}erimental data whereas vibrational- 
rotational close coupling even in tlie simpler e-Ii 2 case poses convergence prob- 
lems of a mopt serious nature, and it still does not include the main effects of 
polarization, which— in a close coupling sense— would require higher lying elec- 
tronic states as well. 

The adiabatic-nuclei amplitudes were calculated from the fixed-nuclei 

amplitudes (the eigenphases of which are those exemplified in Table 1) using 

35 

Nj vibrational functions given by Herman and Wallis. The were 

interpolated, and a fifteen point gaussian quadrature was used to evaluate 

integral (5,3), (Tlie j dependence of the vibrational functions In Ref. 35 was 

surpressed by setting j = 0.) All of these non-IIg amplitudes were combined 

with the amplitudes for the Ilg wave (m = l even parity) from the 

close coupling part of the calculation according the hybrid theory prescription 

36 

above for the particular v and v' given in the results below. 

In Figure 7 we give the integrated vibrationally elastic (v = 0 v' = 0) 

7 

cross section; the experimental result of Golden is shown in the inset. Of 
particular note aside from the substructure itself is the agi^eement of the magni- 
tude of the cross section at the peaks as well as the correct relative size of the 
first two peaks. The remaining peaks are not fully developed and as we have 
indicated this is caused primarily by lack of sufficient C and v coupling due to 
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limitations of machine size. Although the higher pealcs ai’e also small in this 
experiment, they do show up more prominently in other experimental results be- 
low. Tlie theoretical curve does not give any structure below the first pealc at 

11 

1.9eV. Ehrhardt and Willmann have found smooth behaviour here also, so that 
the older experimental structn’'" below 1.8eV would appear to be spurious. The 
major discrepancy with the experiment which is not likely to be altered by further 
coupling is the depth of the calculated cross section between the first two pealts. 
Whether this could be an artifact of the measurement due to lack of sufficient 
energy resolution we cannot say. 

We next turn to differential cross sections involving vibrationally elastic 
scattering. This is the first case we encounter cross term effects in the hybrid- 
theory between resonant (calculated with close coupling) and non resonant (calcu- 
lated by adiabatic-nuclei theory) partial waves. That both cross term and quadratic 
effects are important can be seen by comparing TTigure 8 with Figure 9. Figures 

gives the completely calculated differential cross sections including non- resonant 

11 

waves for three of the energies measured by Ehrhardt and Willmann (which results 
are given in the inset) . In Figure 9 similar results including only the Ilg partial wave 
are given. Not only are the latter spuriously symmetric about 90“ , but the rise 
around 90“ is completely absentboth from the experiment and the full calculation. 

In Figure 10 the full calculation is given for tiie remaining three measured energies. 
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(given in the inset). Note again that the calculation gives absolute values, but 

11 

only the relative values are measured. The shapes appear in satisfactory 
accord with experiment.^^ 

The last comparison we shall give for vibrationally elastic scattering is the 
differential cross section as a function of the energy at various angles (Fig. 11), 
We consider the comparison very satisfactory but the absence of the higher wig- 
gles particularly in the forward directions (due to the lack of sufficient coupling 
in the calculation as explained above) is somewhat more apparent than in the 
integral cross section (Fig. 7). 

We turn next to the inelastic cross sections. In Figures 12 and 13 we give 

the differential values at 0' = 20°, 72° vs. respectively, these being the 

measurements in Reference 11 and Reference 6 respectively. In particular 
6 

Schulz inferred that Ms measurements would reflect the integrated cross sec- 
tions as well. That tMs is so is shown in Figure 14 where the integrated cross 
sections are given and are seen virtually indishinguishable from the 72° curve 

In shape. On the other hand the absolute values wMcb were originally inferred 

0 37 

by Schulz on the basis of summed total cross sections measured by Haas are 

about a factor two smaller than we calculate. Thus the need for direct absolute 

measurements is clear, not only to test the theory but for many atmospheric 

applications, of wMch Reference 1 is one, as well. Although the shape for the 

inelastic v =0 v' = 1 cross section is satisfactory, the detailed agreement 
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rapidly degenerates for higher v\ and again this can be attributed to inadequate 
V and C coupling as noted in going from Figure 4 to Figure 6. We believe how- 
ever that the average magnitudes of the particular cross sections are meaningftil 
for practical applications. 

It should be noted that because of the virtual independence of the nonresonant 

t 

amplitudes on R, essentially all of the inelastic cross sections comes from the 
resonant n wave. This can be clearly seen in the inelastic- differential cross 

o 

sections, given in Figure 16, wMch are quite ssonmetrio about 90°. Although the 
measurements of Ehrliardt and Willmann^^ are not done over sufficiently wide 
angular range to prove the symmetry, the agreement with calculation in terms 
of ratios of forward to minimum to 90° values at the different energies is very 
good. 


Finally we give the momentum transfer cross section. The formula for this 
cross section 

^ (1 -cos0')dJ2' (5.8a) 

is readily integrated from (5.4) to give 


7T 


r2 ^ 
v' 


2 Aq(o, v')“ — Aj(o, v') 
V* L 3 - 


(5.8b) 


where the A(o, v') are given in (5.4b) with the rT„ amplitudes replaced again by 

a 

the close-coupling amplitudes according the prescription (5.7) of the hybrid theory. 
(Note the Aj (o, v') = 0 for v'> 0 for Fig ). The result of the calculation is compared 


25 


J 


with the eiqjeriment in Figure 16. The experiment as is well Icnown does not meas- 
ure directly, but rather infers it by optimizing assumed momentum transfer 
cross sections to fit swarm data as function of applied electric field. As such 
the absence of substructure on the "experimental" result should not be interpreted 
as its absence in fact. (It is virtimlly certain that the substructure must be 
present. ) The experimental curve is seen to envelope the calculated curve in 
tlie resonance region as would be expected. About 25% of the calculated curve 
com.es from, the v' > 0 terms in (5, 8h). However the contribution of the resonances 
does not extend beyond about 4eV, thus the 20% difference in tills energy range 
would appear at this point to remain unexplained. 

Finally it is clear from (5.2) and (5.7) that the hybrid theory can be used to 
calculate simultaneous rotation-vibration excitation. We shall not do that here 
as the programming of the formulae is somewhat more arduous, and there are 
presently no experiments with wliich to compare. (We do intend to perform that 
calculation at a later time.) 

VI. DISCUSSION 

This then completes our hybridization of vibrational close coupling and 
adiabatic-nuclei rotational approximations. In order to complete the a- priori 
theoretical framework for calculation of very narrow (probably Feshbach) 
resonanances, one will have to include rotational coupling as weU. hi that case 
the hybridization will talm place as the second level only: i.e. , tlie rotational 
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and vibrational close-coupling amplitudes for the resonant partial will be merged 
with the adiabatic-nuclei amplitudes for the non-resonantpartial wave . We ex- 
pect to elucidate the formal aspects of this generalization shortly. 

Calculations involving this generalized hybrid theory will certainly be 
arduous. It may well be that for practical purposes other techniques (R-matrix, 
Fredholm determinant, etc.) may be more useful. But to be reliably accurate 
they will have to include the equivalent physics of the hybrid theory. 

The frame-transformation theory^ represents a different mix of the 
above theories. There one ties interior fixed-nuclei calculations to exterior 
close-coupling calculations at a boundary point utilizing appropriate simpli- 
fications for each region. The approach has undoubted utility in molecular photo- 
ionization and electron-molecular ion scattering where the Icnown asymptotic 

41 

Coulomb solutions can be combined with multi-channel quantum defect theory 

to render the resonant structure to be described in terms of a few experimental 
30 

parameters. It should be noted, however, that even here the fixed- and 
adiabatic-nuclei theories can describe the non-resonant structure very 

In the case of scattering fx’om hetero-nuclear molecules (i. e. , those with a 

43 

dipole moment) frame-transformation can be expected to be useful if for no 

other reason than rendering certain cross section finite which would diverge in 

44 

the fixed- and adiabatic-nuclei approximations. However in the case of homo- 
nuclear diatomic molecules the utility of frame transformation is more uncertain, 
because the forces are probably not long-range enough to allow cross section to 
emerge which are suitably insensitive to the matching radius Tq- 
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Table I 

Nonresonant Eigenpbase Sums as a Fimction of R (mod ir ) 


Energy 

(eV) 


Partial 

Wave* 


Interauclear Separation R (in units of Uq) 
1,744393 1.868 2,068 2.268 2,391607 


0.20 

0. 40 
0.60 
0.80 
1.00 
1.20 
1.40 

1. 60 
1.80 
2.00 
2.20 


A 

-0.1921 

B 

-0.0253 

C 

^'.0351 

A 

-0.2960 

B 

-0,0366 

C 

0.0503 

A 

-0.3772 

B 

-0.0505 

C 

0,0572 

A 

-0.4450 

B 

-0.0673 

C 

0.0591 

A 

‘0.5037 

B 

-0.0845 

C 

0.0573 

A 

-0.5555 

B 

-0.1033 

C 

0.0528 

A 

-0.6022 

B 

-0.1229 

C 

0.0466 

A 

-0.6446 

B 

-0.1428 

/*1 

0.0389 

A 

-0.6835 

B 

-0.1630 

C 

0,0303 

A 

-0. 7193 

B 

-0.1833 

C 

0.0210 

A 

-0.7524 

B 

-0.2035 

C 

0.0113 

A 

-0.7833 

B 

-0.2235 

C 

0.0011 


“0.1906 
“0.0645 
0.0295 
-0.2938 
-0.0113 
0.0423 
-0.3747 
0.0211 
0.0474 
-0.4421 
-0.0346 
0. 0479 
-0.5006 
-0.0503 
0.0449 
-0.5524 
-0.0677 
0.0394 
-0.5989 
-0.0862 
0.0323 
-0. 6412 
-0.1054 
0.0239 
-0.6801 
“0.1250 
0.0146 
-0.7158 
-0.1448 
0.0046 
“0.7480 
-0.1647 
-0.0058 
-0.7799 
-0.1845 
-0.0165 


-0.1985 
-0.0037 
0.0288 
-0.3048 
-0.0090 
0.0406 
”0.3878 
- 0. 0202 
0.0448 
-0.4570 
-0.0353 
0.0442 
-0.5170 
-0.0530 
0.0401 
-0.5702 
-0.0723 
0.0337 
-0.6180 
-0.0927 
0.0256 
-0.6616 
-0.1137 
0.0163 
-0.7018 
-0.1352 
0.0062 
-0.7388 
-0.1568 
-0.0045 
-0.7733 
“0.1784 
-0.0155 
-0.8056 
-0.1999 
“0.0269 


-0.1902 
0.0011 
0.0278 
-0.2940 
- 0,0020 
0.0390 
-0.3750 
-0.0115 
0. 0427 
-0.4426 
-0.0253 
0. 0414 
-0.5012 
-0.0419 
0.0368 
-0.5530 
-0.0603 
0.0298 
-0.5997 
-0.0799 
0.0212 
-0.6423 
-0,1003 
0.0114 
-0.6814 
- 0.1211 
0.0009 
-0.7176 
-0.1421 
- 0.0102 
-0.7514 
-0.1632 
-0.0216 
-0.7830 
-0.1843 
-0.0332 


-0.1926 
0.0038 
0.0271 
-0.2972 
0.0012 
0.0379 
-0.3787 
-0,0083 
0.0411 
-0.4467 
-0.0223 
0.0394 
-0.5056 
-0.0393 
0.0343 
-0.5577 
-0.0581 
0.0270 
-0.6047 
-0.0782 
0.0181 
-0.6475 
-0.0990 
0.0081 
-0.6869 
-0.1203 
-0.0027 
-0. 7234 
-0.1417 
-0.0139 
-0.7575 
-0.1632 
-0.0255 
-0.7896 
-0.1846 
-0. 0373 


2.40 
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Table I (Continued) 


Energy 

(eV) 

Partial 

Wave* 

Intemuelear Separation R in Units of ) 

1.744393 

1.868 

2.068 

2.268 

2.391607 

2.60 

A 

-0.8122 

-0.8089 

-0.8359 

-0.8129 

-0.8198 


B 

"0.2434 

-0.2043 

-0.2213 

-0.2051 

-0.2058 


C 

-0.0092 

-0.0274 

-0.0383 

-0.0450 

-0.0491 

2.80 

A 

-0.8393 

-0.8361 

-0.8646 

-0.8409 

-0.8484 


B 

"0.2631 

-0.2239 

-0.2424 

-0.2258 

-0.2268 


C 

-0.0197 

-0.0383 

-0.0498 

-0.0567 

-0.0609 

3,00 

A 

-0.8649 

-0.8617 

-0.8917 

-0.8677 

-0.8756 


B 

-0.2825 

-0.2432 

-0.2633 

-0.2462 

-0.2474 


C 

-0.0302 

-0.0493 

-0.0613 

-0.0684 

-0.0726 

3.20 

A 

-0.8890 

-0.8860 

-0.9175 

-0.8929 

-0.9015 


B 

-0.3016 

-0.2623 

-0.2838 

-0.2664 

-0.2678 


C 

-0.0408 

-0.0603 

-0.0727 

-0.0800 

-0.0841 

3.40 

A 

-0.9119 

-0.9090 

-0.9418 

-0.9174 

-0.9265 


B 

-0.3204 

-0.2812 

-0.3041 

-0.2862 

-0,2878 


C 

-0.0513 

-0.0712 

-0.0839 

-0.0914 

-0.0955 

3.60 

A 

-0.9337 

-0.9307 

-0.9654 

-0.9408 

-0.9505 


B 

-0.3389 

-0.2997 

-0.3239 

-0.3056 

-0.3074 


C 

-0.0617 

-0.0819 

-0.0950 

-0.1026 

-0.1067 

3.80 

A 

-0.9544 

-0.9516 

-0.9880 

-0.9633 

-0.9737 


B 

-0.3570 

-0.3180 

-0.3434 

-0.3247 

“0.3266 


C 

-0.0720 

-0.0925 

-0.1059 

-0.1136 

-0.1177 

4.00 

A 

-0. 9743 

-0.9716 

-1.0096 

-0.9850 

-0.9961 


B 

-0.3749 

-0.3359 

-0.3625 

-0.3434 

-0.3454 


C 

-0.0821 

-0.1030 

-0.1167 

-0.1244 

-0.1284 

4.50 

A 

-1.0203 

-1.0181 

-1.0605 

-1.0363 

-1.0491 


B 

-0.4179 

-0,3792 

-0.4087 

-0.3884 

-0.3907 


C 

-0.1067 

-0.1283 

-0.1425 

-0.1502 

-0. 1541 

5.00 

A 

-1.0620 

-1.0605 

-1.1073 

-1.0838 

-1.0984 


B 

-0.4589 

-0.4204 

-0.4525 

-0.4311 

-0.4332 


C 

-0.1302 

-0.1522 

-0.1609 

-0.1745 

-0.1779 


*Rows A, B, C refer to 2g, eigenphase sums respectively. 
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Figure 1 

Figure 2 

Figure 3 

Figure 4 

Figure 5 
Figure 6 
Figure 7 


Figure 8 
Figure 9 


Figure Captions 

Fixed nuclei IT eigenphase sum (mod tt) for different internuclear 

',7 

separations. Note that k is relative to the ground state energy 

of ^2 which is itself a function R? cf. bottom curve, Fig. 2 . 

Width and position of 11^ resonance vs R. The other results are 

from Birtv;istle and Herzenberg (Ref. 10) and Krauss and Mies 

(Ref. 12). The lowest curve is E« (R) - E ■ 

2 ^2 

Some diagonal and non-diagonal potentials of the vibrational 
close coupling equations vs r. The corresponding fixed-nuclei 
potential is indicated by open circles. 

The Hg contribution to the vibrationally elastic cross section 
Ooo including two and three-partial waves for increasing number 
of vibrational states included in the close coupling expansion - 
Same as Fig. 4 but for (=n]_Q) . 

Same as Fig. 4 but for ^ = *^20^ ' 

The vibrationally elastic scattering in the full hybrid theory. 

In the inset is the experimental results of Golden (Ref. 7) . 

In Ref, 11 the experimental structure below 1.8 eV is not found 
and is considered to be spurious. 

Full hybrid theory calculation of vibrationally elastic differential 
cross sections at three energies. In the inset are experimental 
results of Ehrhardt and Willmann (Ref. 11) . 

Ilg contribution to do^^/dJ^. Note the symmetry and enhancement 
around 90° are not present in the full calculation or the 
experiment; cf. Fig. 8. 
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Figure 10- Full hybrid theory calculation of da^^/dft at remaining three 

energies measured in Ref- 11 which are given in inset. 

2 

Figure 11. da^^^/dfl vs k for various angles. Experimental result of 
Ref. 11 given in inset. 

Figure 12. da^^^/d^^ at 0 = 20° for various excited states. Experimental 

results of Ref. 11 given in inset. 

Figure 13. Same as Fig. 12 for 0 = 72°. In this case the experimental 
result (inset) is that of Schulz, second paper of Ref. 6. 

The ordinate of the inset show the original inferred normal- 
ization; cf. caption of Fig. 14. 

Figure 14. Total vibrational excitation curves- Note similarity of 

o 

shape to 0 = 72 curves Fig, 13 as assumed by Schulz (Ref. 6, 
1964) , from which the experimental result is taken. The 
experimental normalization however is a factor two higher 
than given in Ref. 6 and constitutes a new inferred normal- 
ization (Schulz, 1975 to be published) . In the latter, Schulz 
states that the new normalization may itself be low by a factor 
two, which is in accord with our results. 

Figure 15- Vibrationally inelastic differential cross sections compared 

to experiment (Ref. 11) in inset. Note that the vertical 

, o 

symmetry of the inelastic cross sections around 90 due to 

minuteness of non-Il contribution in these cases, 

y 

Figure 16. Momentum transfer cross section compared to the one inferred 

from swarm experiment by Englehardt, Phelps, and Risk (Ref. 38). 
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